The cx-dependence of transition frequencies for ions Si II, Cr II, Fe II, Ni II, and Zn II 
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We performed accurate calculation of a-dependence (a = e 2 /he) of the transition frequencies for 
ions, which are used in a search for the variation of the fine structure constant a in space-time. 
We use Dirac-Hartree-Fock method as a zero approximation and then the many-body perturbation 
theory and configuration interaction methods to improve the results. An important problem of level 
pseudocrossing (as functions of a) is considered. Near the crossing point the derivative of frequencies 
over a varies strongly (including change of the sign). This makes it very sensitive to the position 
of the crossing point. We proposed a semiempirical solution of the problem which allows to obtain 
accurate results. 
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I. INTRODUCTION 

Recently there was an intensive discussion of the pos- 
sible space-time variation of the fine structure constant 
a — e 2 /he at the cosmological scale. The first evidence 
for such variation has been reported in 0, H H 01 II El 
from the analysis of the astrophysical data. These results 
are to be compared with the number of experimental up- 
per bounds on this variation obtained from other astro- 
physical observations (see, e.g. 0,11 0) and from the 
precision laboratory measurements [TollllL Il2| . Recently 
a number of new laboratory tests have been proposed 
(see, e.g. 13]). The analysis of the microwave back- 
ground radiation can also give some restrictions on time 
variation of a as suggested in 

El El El- Implemen- 
tations of the space-time variation of the fine structure 
constant to the theory of the fundamental interactions 

are discussed e.g. in Refs. Ei ei Ei mm mm see 

also disscussion and references in [3j). 

The most straitforward way to look for the variation of 
a is to measure the ratio of some fine structure interval 
to an optical transition frequency, such as LD{rvp\j% — » 
n P3/2) and uj(n'si/2 — > np 3 / 2 ) |37|]. This ratio can be 
roughly estimated as 0.2 a 2 Z 2 , where Z is the nuclear 
charge j24[. Therefore, any difference in this ratio for 
a laboratory experiment and a measurement for some 
distant astrophysical object can be easily converted into 
the space-time variation of a. However, as it was pointed 
out in 25], one can gain about an order of magnitude 
in the sensitivity to the a-variation by comparing optical 
transitions for different atoms. In this case the frequency 
of each transition can be expanded in a series in a 2 : 



(o) 
= w; 
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(a/a f-l, (lb) 



where cto stands for the laboratory value of the fine struc- 



ture constant. Note, that Eq. (|lafl corresponds to the 
expansion at a = 0, while Eq. Ijlbfl — to the expansion 

(2) 

at a = ao- In both cases parameters and qi appear 
due to relativistic corrections. 

For a fine structure transition the first coefficient on 
the right hand side of l|la|) turns to zero, while for the 
optical transitions it does not. Thus, for the case of a 
fine structure and an optical transition one can write: 



(2) 

Wfe LU f , 4 
— = 3) a 



(2) 



U) 



op 



while for two optical transitions i and k the ratio is: 
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a 2 +0(a 4 ). 



(3) 



(2) 

Quite often the coefficients ui^ for optical transitions are 
about an order of magnitude larger than corresponding 
coefficients for the fine structure transitions cjP (this is 
because the relativistic correction to a ground state elec- 
tron energy is substantially larger than the spin-orbit 
splitting in an excited state [251 l2r|). Therefore, the 
ratio © is, in general, more sensitive to the variation 
of a than the ratio It is also important that the 

signs of coefficients in @ can vary. For example, 
for s-p transitions the relativistic corrections are positive 
while for d-p transitions they are negative. This allows to 
suppress possible systematic errors which "do not know" 
about the signs and magnitude of the relativistic correc- 
tions [25^. On the other hand, for many cases of interest, 
the underlying atomic theory is much more complicated 
for Eq. J2J). In particular, the most difficult case corre- 
sponds to transitions to highly excited states of a multi- 
electron atom, where the spectrum is very dense. And 
this happens to be a typical situation for astrophysical 
spectra, in particular, for large cosmological red shifts. 
Corresponding atomic calculations have to account very 
accurately for the electronic correlations, which may af- 
fect such spectra quite dramatically. 
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TABLE I: Final results for parameters q from Eq. Q for Si II, 
Cr II, Fe II, Ni II, and Zn II. Estimated errors are in brackets. 



Ion 



Sill ' l P° 



Transition 
27 



ojo (cm 



(cm 



/2 



2 Z> 



'ft 



3/2 
/2 



55309.3365 
65500.4492 



520 
50 



Cr II °S B 



/2 



Fell °D. 



^3/2 
6r>o 
r 5/2 
6po 
\7J2_ 



48398 
48491 
48632 



868 
053 
055 



-1360 
-1280 
-1110 



9/2 



Ni II *D, 



U $/2 
^7/2 

*9/2 

^7/2 

^7/2 
6po 
7J2_ 



38458 
38660 
41968 
42114 
42658 
62065 
62171 



9871 

0494 

0642 

8329 

2404 

528 

625 



1330 
1490 
1460 
1590 
1210 
1100 
-1300 



5/2 



Zn II ISi 



b 7/2 

U bl2 
2™ 
^5/2 



57080 
57420 
58493 



373 
013 
071 



-700 
-1400 

-20 



/2 



2t>o 
-<l/2 

f_3/2 



48481 
49355 



077 
002 



1584 
2490 



30) 
30) 



150) 
150) 
150) 



150) 
150) 
150) 
150) 
150) 
300) 
300) 



250) 
250) 
250) 



25) 
25) 



The first calculations of the coefficients q from Eq. 
for the transitions suitable for astronomical and labora- 
tory measurements were done in Refs. [H [2^, [23, l28| . 
Here we present a new and more accurate calculations of 
the coefficients q for the transitions, which are currently 
used in the analysis of the astrophysical data. A full list 
of these transitions was given in [3j. We have not re- 
calculated here the lightest and the most simple atoms 
Mg and Al, for which the previous calculation |25j should 
be sufficiently accurate and focused on more complicated 
ions Si II, Cr II, Fe II, Ni II, and Zn II. Our final results 
for them are given in TablelJ Note, that here we use the 
single parameter q instead of two parameters qi and q<i 
used in the earlier works and q = dtu / dx\ x= Q — qi + 2g 2 - 
Details of the calculations and discussion of the accuracy 
will be given in Sec. IIIII Before that we briefly address 
few theoretical points in Sec. [H] 



II. THEORY 

In order to find parameters q — duj/dx\ x= o in Eq. 
we perform atomic calculations for three values of x: 
X- = —1/8, xq = 0, and x + = 1/8. That allows us 
to determine q: q = 4 (lo(x + ) — lu(x-)) and also estimate 
the second derivative d 2 oj/dx 2 \ x= Q. The large value of 
the latter signals that interaction between levels is strong 
(level pseudocrossing), and there is a risk of large errors. 
For these cases further analysis was done as described 
below. 

a. Relativistic calculations of multi- electron ions. In 
order to accurately account for the dominant relativis- 
tic effects we use the Dirac-Hartree-Fock approximation 
as a starting point for all calculations of atomic spec- 
tra. Though most of the calculations were done for the 



Coulomb potential, we have also estimated Breit correc- 
tions by including the magnetic par t of the Breit inter- 
action in the self-consistent field |29|. 

The ions we are dealing with in this paper have from 
one to nine electrons in the open shells. For one valence 
electron in Zn II the Dirac-Fock V N ~ X approximation 
already gives rather good results. On the next step the 
core-valence correlations can be accounted for by means 
of the many-body perturbation theory (MBPT). Already 
the second order MBPT correction allows to reproduce 
the spectrum with the accuracy, better than 1%, which 
is more than sufficient for our current purposes. 

Other ions of interest to us have at least three valence 
electrons. Here the dominant correlation correction to 
transition frequencies corresponds to the valence- valence 
correlations. This type of correlations can be accounted 
for with configuration interaction (CI) method. If neces- 
sary, the core- valence correlations can be included within 
combined CI+MBPT technique [13 ■ The latter usually 
provides an accuracy of the order of 1% or better for the 
lower part of the spectra of atoms and ions with two or 
three valence electrons 32]. However, the accu- 

racy of ab initio methods decreases with the number of 
valence electrons and with excitation energy. Indeed, for 
a large number of valence electrons and/or sufficiently 
high excitation energy the spectrum becomes dense and 
the levels with the same exact quantum numbers strongly 
interact with each other. The part of the spectrum of Fe 
II above 55000 cm -1 and, to a somewhat lesser extent, 
the spectrum of Ni II represent this situation. There- 
fore, for these ions we developed a semiempirical fitting 
procedure, which is described below. 

In order to have additional control of the accuracy of 
our CI we performed calculations for most of the ions 
with two different comput er p ackages. One package was 
used earlier in Refs. l3Ct l32t |33 | and another one was 
used in Refs. 

U El S imiinEni The former pack- 
age allows to construct flexible basis sets and optimize 
configuration space, while the latter allows for a larger 
CI space as it works with the block of the Hamiltonian 
matrix, which corresponds to a particular total angular 
momentum of atom J. When there were no significant 
difference between two calculations, we only give results 
obtained with the first package. Nevertheless, our final 
results presented in Table [U are based on both calcula- 
tions. 

b. Semiempirical treatment of the strong interaction 
of levels: pseudo-crossing. In the nonrelativistic limit 
a — > 0, all multi-electron states are accurately described 
by the LS-coupling scheme: E a ^o = -E Pl n,L,s,,/, where 
p = ±1 is the parity and n numerates levels with the 
same p, L,S, and J. For sufficiently small values of a the 
LS-coupling holds and the energy has the form: 



Ep, n ,L,S,J - E p,i,L,S + y — J ^Cp,™,L,S (4) 

+\a p ,u,l,s [J(J + 1) - L(L + 1) - S{S + 1)A , 
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where the first term in the parentheses gives the slope 
for the centre of the multiplct and the second term gives 
the fine structure. With growing a the multiplets start 
to overlap and when the levels with the same p and J 
come close, the pseudo-crossing takes place. 

Near the pseudo-crossing the slope of the energy curves 
changes dramatically. If such crossing takes place at x « 
0, where x is defined by Eq. i.e. near the physical 
value of a, it can cause significant uncertainty in the 
values of parameters q. 

Let us first analyze the behaviour of the slopes q(x) 
in the vicinity of the pseudo-crossing in the two-level ap- 
proximation. Consider two levels E\ and E 2 which cross 

0jt X — X f • 



Ei = qx(x — x c ), 
E 2 = q 2 (x-x c ). 



(5a) 
(5b) 



If the interaction matrix element between these two levels 
is V, the exact adiabatic levels will be 



E a ,b = \\Ml +12){X-X C ) 

± V(qi-q2) 2 (x- Xc r + w 2 ). (6) 



It is easy now to calculate the energy derivative in respect 
to x in terms of the mixing angle (f> between unperturbed 
states I and 2: 



dE a , 
dx 



?1,2 



52,1- 



(7) 



Note, that at the crossing the angle <f> varies from 
on one side through 7r/4 in the centre to 7r/2 on 
the other side, which leads to the change of the slope 
q a (x) = dE a /dx from qi through (qi + q 2 )/2 to q 2 . The 
narrow crossings with small V are particularly danger- 
ous, as the slopes change very rapidly within the interval 
Ax ~ V/\qi — q 2 \. Then, even small errors in the posi- 
tion of the crossing point x Cl or the value of V can cause 
large errors in q a ,b- In this model we assume that non- 
diagonal term V = const. For the real atom V cx a 2 . 
However, if the crossing region Ax < 1, we can neglect 
the dependence of 7 on a. 

c. Semiempirical treatment of the strong interaction 
of levels: multi-level case. Eq. (Q can be easily gener- 
alized to a multi-level case as it simply gives the slope 
of a physical level a as a weighted average of the mixed 
levels. Thus, if the level a can be expressed as a linear 
combination of some unperturbed LS'-states ipL n ,s n '- 



\a) = Y,C n \ip Lnt s n }, 



the resultant slope q a is given by: 

9a = Yl C 2 q n 



(8) 



(9) 



Here again we neglect weak dependence of interaction V 
on x in comparison to strong dependence of C 2 on x near 
crossing points. 

Eq. J!|J) allows to improve ab initio coefficients q if we 
can find the expansion coefficients C n in Eq. JSJ). That 
can be done, for example, by fitting (/-factors. The mag- 
netic moment operator /x = go(L + 2S) is diagonal in L 
and S and, therefore, does not mix different LS'-states. 
Thus, in the LS'-basis the resultant g-factor for the state 
a has exactly the same form as q a : 



9a 



(10) 



If the experimental ^-factors are known, one can use 
Eq. (IKjp to find weights C 2 and, then find the corrected 
values of the slopes q a . 

Sometimes, the experimental data on g-factors are in- 
complete. Than, one can still use a simplified version of 
Eqs. © and CHI): 



g a = C 2 g° + (1-C 2 )g, =► C 2 = 9 -^4, 
9a = C 2 q a a + (1-C 2 )q. 



(Ha) 
(Hb) 



C 2 here is the weight of the dominant LS-\eve\ in the ex- 
perimental one, and the bar means the averaging over the 
admixing levels. Of course, there is some arbitrariness in 
calculation of averages g and q. However, the advantage 
of Eqs. (f 1 1 1> is that only one experimental g-factor is re- 
quired. 



III. DETAILS OF THE CALCULATION AND 
RESULTS 

As we mentioned above, we performed calculations of 
energy levels for three values of the parameter x: X- = 
— 1/8, x = 0, and x + = 1/8. All three calculations were 
done at exactly same level of approximation, to minimize 
the error caused by the incompleteness of the basis sets 
and configuration sets. From these calculations we found 
two approximations for q: g_ = 8(u(xo) — uj(x^)) and 
g+ = 8(lj(x+) —uj(xq)). If there were problems with level 
identification we performed additional calculation for x = 
0.01, where the LS-coupling should be very accurate and 
identification is straitforward. The noticeable difference 
between q_ and q + signaled the possibility of the level 
crossing. In these cases we applied the semiempirical 
procedure described in Sec.|H]to find the corrected values 
for q; otherwise, we simply took the average: q = {q+ + 
q-)/2. 



Zn II 



Zn II has the ground state configuration [Is . . . 3d ] 4s 
and we are interested in the 4s — > Apj transitions. As 
the theory here is much simpler than for other ions, we 
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TABLE II: Transition frequencies and parameters q for Zn II 
(in cm -1 ). Calculations were done in four different ap- 
proximations: Dirac-Hartree-Fock-Coulomb (DHFC), Dirac- 
Hartree-Fock-Coulomb-Breit (DHFCB), Brueckner-Coulomb 
(BC), and Brueckner-Coulomb-Breit (BCB). 



Transition 


Exper. DHFC 


DHFCB 


BC 


BCB 


4s 1/2 -> 4pi/ 2 
-> 4p 3/2 


transition frequencies 
48481.077 44610.1 44608.1 
49355.002 45346.9 45330.0 


48391.2 
49263.8 


48389.4 
49244.6 


4s 1/2 -> 4p 1/2 
-» 4p 3/2 


parameters q = 
1362 
2129 


(q+ + q- 

1359 
2109 


)/2 
1594 
2500 


1590 
2479 



TABLE III: Transition frequencies ui from the ground state 
2 P°/ 2 , fine structure splitting Afs, and parameters q± for Si II 
(in cm -1 ). 



r 3/2 
Pi/ 2 
"3/2 
P 5/2 



2 '5'l/2 



3/2 
5/2 



Experiment [35] 




Theory 








Afs 


LU 


Afs 


q- 


q+ 


287 


287 


293 


293 


295 


291 


44080 




41643 




453 


451 


44191 


111 


41754 


111 


565 


564 


44364 


174 


41935 


181 


746 


744 


55304 




54655 




509 


507 


55320 


16 


54675 


20 


530 


530 


65495 




65148 




40 


39 



used Zn II to study the importance of the core-valence 
correlation correction and Breit correction to the slopes 
q. The former correction was calculated in Brueckner 
approximation: 

(flDHF + £(#)) * = E *> ( 12 ) 

with the self-energy operator T,(E) calculated in the sec- 
ond order of MBPT (the perturbation here is the differ- 
ence between the exact and Dirac-Hartree-Fock Hamil- 
tonians, V = H — Huhf)- The i?DHF was calculated 
with the magnetic part of the Breit operator included 
self-consistently. The retardation part of the Breit op- 
erator is known to be significantly smaller p9l | and we 
completely neglected it here. 

The results of our calculations of the frequencies uo and 
the slopes q for two transitions 4s — ► 4pj, j = 1/2,3/2 
are given in Table ITU One can see, that both Brueckner- 
Coulomb and Brueckner-Coulomb-Breit approximations 
give very good transition frequencies, accurate to 0.2%, 
though the latter slightly underestimates the fine split- 
ting. Breit correction to the parameters q does not exceed 
1%, while core-valence correlations account for the 17% 
correction. 

In Table ITTI we did not give separately the values of q±. 
The difference between them is close to 1%. Indeed, in 
the absence of close interacting levels the dependence of 
q on x arise from the corrections to the energy of the 
order of a 4 Z 4 , which are very small. 

B. Si II 

Si II has three valence electrons and the ground state 
configuration [Is 2 . . . 2p 6 ]3s 2 3p. Excited configurations 
of interest are 3s3p 2 and 3s 2 4s. We made the CI cal- 
culation in the Coulomb approximation on the basis set, 
which included Is — 8s, 2p — 8p, 3d — 8d, and 4/, 5/ 
orbitals, which we denote as the basis set [8spd5f]. Note, 
that we use virtual orbitals, which are localized within 
the atom |34j], rather than Dirac-Fock ones. This pro- 
vides fast convergence. CI included all single-double (SD) 
and partly triple excitations from three valence configu- 
rations listed above. The results of these calculations are 
given in Tabic ITTT1 



Like in Zn, the left and write derivatives q- and q+ 
are close to each other, and all levels with equal exact 
quantum numbers are well separated. The astrophysical 
data exist for the levels 2 Si/ 2 and 2 D 5 / 2 - The former cor- 
responds to the 3p — > 4s transition and has small slope 
g, while the latter corresponds to the 3s — > 3p transition 
and has much larger positive q. That is in agreement with 
the fact, that relativistic corrections to the energy usu- 
ally decrease with the principle quantum number n and 
with the orbital quantum number I. Therefore, for the 
ns — > np transition one should expect large and positive 
g, while for np — > [n + l)s, there should be large can- 
cellation of relativistic corrections to upper and to lower 
levels, resulting in smaller q (see disscussion in poll^ ). 
The dominant correction to our results should be from 
the core- valence correlations. In the recent calculations of 
Mg, which has the same core as Si II, the core- valence cor- 
rections to transition frequencies were found to be about 
4% [33l ■ We conservatively estimate corresponding 
correction to q to be 6% of the larger q, i.e. 30 cm -1 . 



C. Cr II 

Cr II has the ground state configuration [Is 2 . . . 3p 6 ]3d 5 
with five valence electrons. The astrophysical data cor- 
respond to the 3d — > Ap transition, for which one may 
expect negative value of q. CI calculations here are much 
more complicated, than for Si II. There is strong relax- 
ation of the 3d shell in the discussed transition, which 
requires more basic d-orbitals. Therefore, we used the 
[6sp9d6f] basis set. In CI we included only single and 
double (SD) excitations. Some of the triple, quadruple, 
and octuple excitations were accounted for by means of 
the second order perturbation theory. It was found that 
corresponding corrections to transition frequencies were 
of the order of few percent, and were even smaller for pa- 
rameters q. In general, these corrections did not improve 
the agreement with the experiment, so we present only 
CI results in Table ITVI 

As we mentioned above, there is strong relaxation of 
the 3d-shell in the 3d —> Ap transition. We were not 
able to saturate CI space and completely account for this 
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TABLE IV: Transition frequencies uj from the ground state 
6 S 5 /2, fine structure splitting Aps, and parameters q for Cr II 
(in cm ). CI single-double approximation was used for the 
Coulomb-Breit interaction. 



Experiment 

LU Aps 



Theory 
A PS 



q+ 



»D 



/2 



7/2 
) 

/2 



^3/2 

^5/2 
61710 
^7/2 

-^9/2 
6/710 

-^11/2 
6po 

■M/2 
6po 

■ f 2/2 
6po 

r 3/2 



12148 

12304 
46824 
46906 
47041 
47228 
47465 
47752 
48399 
48491 
48632 



13123 
156 13289 

47163 
82 47244 
135 47378 
187 47565 
237 47803 
287 48091 

48684 
92 48790 
141 48947 



165 

81 
134 
187 

238 
288 

106 
157 



-2314 
-2153 
-1798 
-1715 
-1579 
-1387 
-1148 
-862 
-1364 
-1278 
-1108 



relaxation. Because of that, we estimate the error for q 
here to be close to 10%. 

We have seen before for Zn II and Si II, that in the 
absence of level-crossing the difference between q + and 
g_ is smaller than other theoretical uncertainties. In Cr 
II there are no close levels which may interact with each 
other, so in the calculation presented in Table HVI we de- 
termined only the right derivative q+. In calculations 
with different basis sets we checked that the difference 
between q + and g_ is much smaller than the given above 
theoretical error (see Table UJ. 



D. Fe II 

Fe II ion has 7 valence electrons in configuration 3d 6 As 
and represents the most complicated case. The astro- 
physical data includes 5 lines in the band 38000 cm _1 - 
43000 cm -1 and two lines with the frequency close to 
62000 cm -1 . The first band consists of three close, but 
separated multiplets with a regular fine structure split- 
tings. The 62000 cm _1 band is completely different as the 
multiplets here strongly overlap and fine structure inter- 
vals are irregular Characteristic distance between 
the levels with identical exact quantum numbers is few 
hundred cm -1 , which is comparable to the fine structure 
splittings. This means that the levels strongly interact 
and even their identification may be a problem. 

In fact, in Moore Tables [3f| one of the multiplets of 
interest, namely y e P°, is erroneously assign to the con- 
figuration 3d 6 { 7 S)4p. It is an obvious misprint, as there 
is no term S for configuration 3d 6 . This term appears, 
however, in the configuration 3c? 5 and the correct assign- 
ment of this multiplet should be 3d 5 ( 7 S')4s4p. This as- 
signment is in agreement with our calculations and with 
the experimental g-factor of the level with J = 7/2. We 
checked that all close levels of the configuration 3d 6 Ap 
have significantly smaller g-factors. 

This reassignment has dramatic consequences in terms 




FIG. 1: Examples of typical interaction of levels in the upper 
band of Fe II. Levels are shown in arbitrary units as function 
of (a/ao) 2 = x + 1. Levels of configuration 3d 6 4p have similar 
slopes and strongly interact with each other. That causes 
wide pseudo-crossings, similar to one shown on the left side 
of the plot. The level 6 Py 2 of the configuration 3d 5 4s4p moves 
in the opposite direction. A series of sharp pseudo-crossings 
takes place near the physical value of a, marked by a vertical 
dotted line. 



of the corresponding parameter q as configurations 3<i 6 4p 
(4s — 4p transition from the ground state) and 3d 5 As4p 
(3d — 4p transition) move in the opposite directions from 
the ground state configuration 3d 6 4s when x is changed. 
It also causes a number of pseudo-crossings to occur right 
in the vicinity of x = (see Fig.^l. 

CI calculations for Fe II were done on the basis set 
[6spdf] in the SD approximation (see Table W\ . Triple 
excitations were included within second order perturba- 
tion theory and corresponding corrections were found to 
be relatively small. One can see from Table[3that for the 
lower band both frequencies and g-factors are reproduced 
rather accurately. 

The first anomaly takes place at 44000 cm -1 , where 
the levels A Dy 2 and A Fy 2 appear in the reverse order. 
Theoretical g-factors are also much further from LS val- 
ues (1.429 and 1.238). That means that theoretical levels 
are at pseudo-crossing, while experimental levels already 
passed it. Indeed, calculations for x — 1/8 show that the 
right order of levels is restored, though the (/-factors are 
still too far from LS values. 

The second anomaly corresponds to the band above 
60000 cm -1 . Here the order of calculated levels differs 
from that of the experimental ones. Note, that for this 
band only levels of negative parity with J = 7/2 are 
given in Table [3 Thus, all of them can interact with 
each other. Let us estimate, how this interaction can 
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TABLE V: Transition frequencies uj from the ground state 
6 Dg/ 2 , ^-factors, and parameters q± for Fe II (in cm -1 ). 



Experiment 



Theory 



^9/2 
6rjo 
LJ 7/2 
67710 
•^11/2 
6770 
^9/2 
67-10 
t 7/1 
&JDO 



1 7 



4 D 

4 F 7 C 



/2 
7/2 



/2 
Spo 
^7/2 

^7/2 
4rjo 
n 7/2 

47710 

t 7/2 

> 

/2 



6po 
^7 

2^0 
^7/2 



a; 


9 


UJ 


9 




9- 


1+ 


oo4oy 


1.04Z 






1.000 


looy 


I0O0 


OoOOU 


1.054 


Oo004 


1.000 


l.Oo ( 


lozz 


1 cin 
101U 


41yoo 




/I 1 Qfi -1 

41oo4 
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affect the slopes q. 

Five levels from this band belong to configuration 
3d 6 Ap and have close slopes with the average q = 
1360 cm -1 . Only the level 4 Fy 2 has the slope, which 

is 300 cm -1 smaller, than the average. The remaining 
level 6 Py 2 belongs to configuration 3<i 5 4s4p and has the 
slope of the opposite sign qi = —1519 cm -1 . Its abso- 
lute value is 500 cm -1 smaller, than for the level 8 Py 2 

of the same configuration 3d 5 4sAp. That suggests that 
the levels 4 F°, 2 and 6 Py 2 strongly interact with each 
other. This is also in agreement with the fact, that 
these levels are the closest neighbors both experimen- 
tally and theoretically and that they cross somewhere 
between x- and x. There is also strong interaction be- 
tween the levels ^Gy/ 2 ) Fy 2 , an< ^ ^7/2- That can be 
seen if one calculates the scalar products (overlaps) be- 
tween corresponding wave functions for different values 
of x, such as: (i(x~)\k(x+)}. For weekly interacting lev- 
els (i(x-)\k(xo)) pa (i(x-)\k(x + )) w Si t k, so large non- 
diagonal matrix elements signal, that corresponding lev- 
els interact. 

Interaction of levels \iy 2 , 4 Fy 2 i anc ^ ^7/2 does not 
affect the slopes q as strongly, as the interaction of 4 F. y 2 
and e Py 2 , so we can account for the former in a less 
accurate way, but it is important to include the latter as 
accurately as possible. 

The level Py 2 interacts with some linear combination 
of levels 2 Gy 21 4 Fy 2 , and 4 Dy 2 . The slopes and g-factors 
of the latter are relatively close to each other, so we can 
simply take the average for all three: 

3 = 1.185; g = 1297. (13) 

Now we can use experimental g-factor of the state 6 P£ , 2 



FIG. 2: Dependence of the odd levels of Ni II on (a/ao) 2 = 
x + 1. Solid lines correspond to J — 5/2 and dashed lines to 
J = 7/2. The experimental positions of the lines are shown as 
short horizontal lines and are all shifted by 1000 cm" 1 . The 
order of levels from bottom up: 4 -D? /2i5/2 , 4 G?/ 2 , 5 /2. 4 ^7°/2,5/2. 

C?/2! -P5/2! -^5/2) an d F 5 j 2 . 

and Eq. (|11|) to determine the mixing: 

C 2 = ^-l = 0.937, (14) 
1.713 - g 

q( 6 P y2 ) = -1342. (15) 

Eq. I|15|) corresponds to the correction Sq = +177. There- 
fore, for the closest level 4 F°, 2 this model gives an esti- 
mate: 

q(%° /2 ) = q - 8q = 1120. (16) 

Eqs. I|15|l and (|16fl show that correction for the mix- 
ing is not very large. That corresponds to the fact that 
experimental g-factor of the level 6 Py 2 is significantly 
larger than any g-factors of the levels of the configuration 
3d 6 Ap. Thus, the interaction for this level is relatively 
small. On the contrary, the levels of the configuration 
3d 6 Ap strongly interact with each other, but correspond- 
ing changes of the slopes are also relatively small (since 
the q values for these strongly interacting levels are ap- 
proximately the same). 

We estimate the accuracy of our calculations for the 
lower band of Fe II to be about 150 cm , and approxi- 
mately 300 cm -1 for the values l|15|) and l|16|l . 

E. Ni II 

Ni II has the ground state configuration 3d 9 . The spec- 
trum is somewhat simpler, than for Fe II. There are also 
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TABLE VI: Transition frequencies uj from the ground state 
2 D 5 /2, g~ factors, and parameters q± for Ni II (in cm -1 ). 



4 n° 

4 n° 

X "7/2 
4f~io 
"5/2 

4 7?° 

t 7/2 
47710 

^5/2 
2s-to 

"7/2 

^7/2 

2 n° 
2 cm 



Experiment 



Theory 

9 g(LS) 



q+ 



1507 
51558 
52739 
54263 
55019 
55418 
56075 
56372 
57080 
57420 
58493 



1.420 

1.356 

1.02 

0.616 

1.184 

0.985 

0.940 

1.154 

1.116 

0.946 



1579 
50415 
51640 
53150 
53953 
54323 
55063 
55284 
56067 
56520 
57589 



1.423 
1.360 
1.016 
0.617 
1.183 
0.986 
0.933 
1.128 
1.108 
0.959 



0.800 
1.429 
1.371 
0.984 
0.571 
1.238 
1.029 
0.889 
1.143 
1.200 
0.857 



1559 
-2405 
-1217 
-1334 
-370 
-1104 
-332 

-60 
-911 
-1419 

-35 



1552 
-2425 
-1245 
-1387 
-418 
-1124 
-334 
-188 
-713 
-1438 
-5 



pseudo-crossings here, but they either lie far from x = 0, 
or are rather wide. That makes their treatment slightly 
easier. Nevertheless, our results significantly differ from 
previous calculations [2^ ■ 

CI calculations were done for the Coulomb potential 
and included SD and partly triple excitations on the ba- 
sis set [5spdf]. We calculated 5 lower odd levels with 
J = 5/2 and 5 with J = 7/2 for X-, xq, and x+, and used 
parabolic extrapolation for the interval —0.4 < x < +0.3 
(see Fig. It is seen that the theory accurately repro- 
duce relative positions of all levels. An overall agreement 
between the theory and the experiment becomes close to 
perfect if all experimental levels are shifted by 1000 cm -1 
down, as it is done in Fig. El Note, that this shift consti- 
tutes only 2% of the average transition frequency. 

Calculated ^-factors are generally in agreement with 
the experiment |35j and noticeably different from the 
pure LS- values (see Table IVT|) . However, for the level 
F°, 2 theoretical <?-factor is smaller than the LS value, 
while experimental one is larger. There are no nearby 
levels who may mix to this one and move (/-factor closer 
to experiment. On the other hand, the difference with 
experiment is only 2% and may be within experimental 
accuracy. 

Fig. |21 shows that the levels ^7/2 an d 2 F°, 2 cross at 
x w 0.3 and they already strongly interact at x = 0. 
Theoretical splitting for these levels is 10% larger than 
experimental one. Thus, they are in fact even closer to 
the crossing point than is predicted by the theory. The 
experimental splitting is equal to the theoretical one for 
larger value of a corresponding to x ~ 0.15. At x = 0.15 
the slopes of these levels are —265 and —590, and for x 
they are —124 and —812 correspondingly. Note, that 
the sum of the slopes at x — 0.15 differs by 80 cm -1 from 
the sum at x = 0. According to Eq. J7J for a two-level 
system the sum is constant. This means that these two 



levels repel from the lower lying level A Fy 2 . Taking this 
analysis into account we suggest an average between x = 
and x = 0.15 as our final value: q( 2 F° /2 ) = -700(250). 



Conclusions 

In this paper we present new refined calculations of 
the parameters q, which determine a-dependence of the 
transition frequencies for a number of ions used in the 
astrophysical search for a-variation. These ions appear 
to be very different from the theoretical point of view. 
Because of that we had to use different methods and dif- 
ferent levels of approximation for them. The final accu- 
racy of our results differs not only for different ions, but 
also for different transitions. 

The simplest system is Zn II, which has one valence 
electron. On the other hand, this is the heaviest ion and 
it has the largest core, which includes 3ci 10 -shell. That 
gave us the opportunity to study corrections to q from the 
core-valence correlations and from Breit interaction. We 
found the former to be about 17% and the latter to be less 
than 1%. For lighter ions Breit interaction should be even 
smaller and can be safely neglected. Other ions also have 
much smaller and more rigid cores, so one might expect 
that core- valence correlations are few times weaker there 
in comparison to Zn. That allows us to neglect core- 
valence correlations for all other ions discussed in this 
paper. 

Si II has the smallest core Is 2 . . . 2p 6 and three valence 
electrons. For neutral Mg, which has the same core, the 
core-valence corrections to the 3s — » 3p transition fre- 
quencies were found to be about 4% |33, Hi| . CI calcula- 
tion for Si II is relatively simple and the errors associated 
with incompleteness of CI space are small. Thus, our es- 
timate of the accuracy for Si on 6% level seems to be 
rather conservative. 

Cr, Fe, and Ni have the core Is 2 . . . 3p 6 and the core 
excitation energy varies from 2 a.u. for Cr II to 2.6 a.u. 
for Ni II. In comparison, the core excitation energy for 
Zn II is 0.9 a.u. Therefore, we estimate the core-valence 
correlation corrections for these ions to be at least two 
times smaller, than for Zn II. 

Additional error here is associated with incomplete- 
ness of the CI space. These ions have from 5 to 9 valence 
electrons and CI space can not be saturated. To estimate 
corresponding uncertainty we performed several calcula- 
tions for each ion using different basis sets and two dif- 
ferent computer packages described in Sec.[n] The basic 
Dirac-Hartree-Fock orbitals were calculated for different 
configurations (for example, for the ground state config- 
uration and for excited state configuration, etc.). 

Supplementary information on the accuracy of our cal- 
culations can be obtained from comparison of calculated 
spectra and <?-factors with experimental values. The later 
appear to be very important as they give information 
about electron coupling, which depends on relativistic 
corrections and on interaction between LS'-multiplets. 
Our results for Cr II appear to be very close for dif- 
ferent calculations and are in good agreement with the 
experiment both in terms of the gross level structure and 
spin-orbit splittings (see Table llVjl . so we estimate our 
final error here to be about 10 - 12%. 
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The largest theoretical uncertainties appear for Fe II 
and Ni II where the number of valence electrons is largest 
and the interaction of levels is strongest. Here we had 
to include semi-empirical fits to improve the agreement 
between the theory and the experiment. We took into 
account the size of these semi-empirical corrections in 
estimates of the accuracy of the calculated values of q. 

The final results are presented in Tabled Note again, 
that they are based on several independent calculations 
performed using two different computer codes. Some of 



the intermediate results are given in Table ITT1- Table IVTl 



Acknowledgments 

This work is supported by Australian Research Coun- 
cil. One of us (MK) thanks UNSW for hospitality and 
acknowledges support from the Gordon Godfrey Fund. 



[1] J. K. Webb, V. V. Flambaum, C. W. Churchill, M. J. 

Drinkwater, and J. D. Barrow, Physical Review Letters 

82, 884 (1999). 
[2] J. K. Webb, M. T. Murphy, V. V. Flambaum, V. A. 

Dzuba, J. D. Barrow, C. W. Churchill, J. X. Prochaska, 

and A. M. Wolfe, Phys. Rev. Lett. 87, 091301 (2001). 
[3] M. T. Murphy, J. K. Webb, V. V. Flambaum, V. A. 

Dzuba, C. W. Churchill, J. X. Prochaska, J. D. Barrow, 

and A. M. Wolfe, Mon. Not. R. Astron. Soc. 327, 1208 

(2001). 

[4] M. T. Murphy, J. K. Webb, V. V. Flambaum, C. W. 

Churchill, and J. X. Prochaska, Mon. Not. R. Astron. 

Soc. 327, 1237 (2001). 
[5] M. T. Murphy, J. K. Webb, V. V. Flambaum, M. J. 

Drinkwater, F. Combes, and T. Wiklind, Mon. Not. R. 

Astron. Soc. 327, 1244 (2001). 
[6] M. T. Murphy, J. K. Webb, V. V. Flambaum, C. W. 

Churchill, and J. X. Prochaska, Mon. Not. R. Astron. 

Soc. 327, 1223 (2001). 
[7] F. Combes and T. Wiklind, Astron. Soc. Pacific. Astro- 
nomical Society of the Pacific Conference Series 156, 210 

(1999). 

[8] A. V. Ivanchik, A. Y. Potekhin, and D. A. Varshalovich, 
Astronomy & Astrophysics 343, 439 (1999). 

[9] D. A. Varshalovich, A. V. Ivanchik, and A. Y. Potekhin, 
Zhurnal Tekhnicheskoi Fiziki 69, 1 (1999). 
[10] J. D. Prestage, R. L. Tjoelker, and L. Maleki, AIP. Amer- 
ican Institute of Physics Conference Proceedings (457), 
357 (1999). 

[11] Y. Sortais, S. Bize, C. Nicolas, M. Santos, C. Mandache, 
G. Santarelli, C. Salomon, and A. Clairon, Proceedings 
of the 1999 Joint Meeting of the European Frequency 
and Time Forum and the IEEE International Frequency 
Control Symposium. IEEE. 1, 34 (1999). 

[12] Y. Sortais, S. Bize, C. Nicolas, G. Santarelli, G. S. C. Sa- 
lomon, and A. Clairon, IEEE Transactions on Ultrasonics 
Ferroelectrics & Frequency Control 47, 1093 (2000). 

[13] C. Braxmaier, O. Pradl, H. Muller, A. Peters, J. Mlynek, 
V. Loriette, and S. Schiller, Physical Review D 64, 
042001 (2001). 

[14] M. Kaplinghat, R. J. Scherrer, and M. S. Turner, Physi- 
cal Review D 60, 023516 (1999). 
[15] S. Hannestad, Physical Review D 60, 023515 (1999). 
[16] J. Kujat and R. J. Scherrer, Physical Review D 62, 



023510 (2000). 

[17] C. Alvarez and R. B. Mann, Physical Review D 55, 1732 

(1997) . 

[18] L. Bergstrom, S. Iguri, and H. Rubinstein, Physical Re- 
view D 60, 045005 (1999). 

[19] R. W. Kuhne, Modern Physics Letters A 14, 1917 (1999). 

[20] N. Chamoun, S. J. Landau, and H. Vucetich, Physics 
Letters B 504, 1 (2001). 

[21] D. Kabat and A. Rajaraman, Physics Letters B 516, 383 
(2001). 

[22] Y. Fujii, Phys. Rev. D 62, 064004 (2000). 

[23] H. Sandvik, J. D. Barrow, and J. Magueijo, A 

simple varying-alpha cosmology (2001), e-print: 

|asTro-ph/0107512l 
[24] I. I. Sobelman, Atomic spectra and radiative transitions 

(Springer- Verlag, Berlin, 1979). 
[25] V. A. Dzuba, V. V. Flambaum, and J. K. Webb, Physical 

Review Letters 82, 888 (1999). 
[26] V. A. Dzuba, V. V. Flambaum, and J. K. Webb, Physical 

Review A 59, 230 (1999). 
[27] V. A. Dzuba and V. V. Flambaum, Physical Review A 

61, 034502 (2000). 
[28] V. A. Dzuba, V. V. Flambaum, M. T. Murphy, and J. K. 

Webb, Phys. Rev. A 63, 042509 (2001). 
[29] E. Lindroth, A.-M. Martensson-Pendrill, A. Ynnerman, 

and P. Oster, J. Phys. B 22, 2447 (1989). 
[30] V. A. Dzuba, V. V. Flambaum, and M. G. Kozlov, Phys. 

Rev. A 54, 3948 (1996). 
[31] V. A. Dzuba and W. R. Johnson, Phys. Rev. A 57, 2459 

(1998) . 

[32] M. G. Kozlov and S. G. Porsev, Opt. Spectrosc. 87, 352 

(1999) . 

[33] S. G. Porsev, M. G. Kozlov, Y. G. Rakhlina, and A. Dere- 

vianko, Phys. Rev. A 64, 012508 (2001). 
[34] P. Bogdanovich and G. Zukauskas, Sov. Phys. Collection 

23(5), 18 (1983). 
[35] C. E. Moore, Atomic Energy Levels, vol. 1-3 (Natl. Bur. 

Stand. (US), Circ. No. 467, Washington, 1958). 
[36] S. G. Porsev, M. G. Kozlov, and Y. G. Rahlina, JETP 

Lett. 72, 595 (2000). 
[37] In fact, the frequency uj(np 1 / 2 — * np 3 / 2 ) is not measured 

directly, but is found as a difference: ui(n's 1 / 2 — * rip :i f 2 ) — 

uj(n's 1/2 -> np 1/2 ). 



